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Exclusive electroproduction of hard dijets can be described within the collinear factorization. 
This process has clear experimental signature and provides one with an interesting alternative 
venue to test QCD description of hard diffractive processes and extract information on generalized 
nucleon parton distributions. In this work we present detailed leading-order QCD calculations of 
the relevant cross sections, including longitudinal momentum fraction distribution of the dijets and 
their azimuthal angle dependence. 
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I. INTRODUCTION 

The QCD description of hard diffraction presents 
an interesting challenge at the crossroads of soft 
and hard physics and appeals for a synthesis of var- 
ious theoretical approaches. In particular diffrac- 
tive exclusive dijet production in deep-inelastic 
lepton-nucleon scattering has attracted consider- 
able attention |3 EJ El This process can in- 
tuitively be visualized as the incident virtual pho- 
ton disintegration into a quark-antiquark pair with 
large and opposite transverse momenta 



e(l)p(p) -> e(l')q( qi )q(q 2 ) P ( P ') 



(1) 



Here I, I' and p, p' are the initial and the final 
momenta of the lepton and the nucleon, respec- 
tively, while q\ and qi are the jet momenta which 
to the tree level accuracy can be identified with 
the momenta of the outgoing quark and antiquark, 
see Fig. ^ In the following discussion we will use 
conventional variables 



= 1-1', 

Q 2 

2(p ■ q) ' 



y 



(q-p) 
(i-p) 



(2) 



and neglect the proton mass whenever it is possible. 

Our interest to the reaction in (JTJ is twofold. 
First, this process provides one with a sensitive 
probe of the generalized gluon parton distribution 
(GPD) in the nucleon @, 0, 

see e.g. [8j. In 
this quality, the dijet production with large invari- 
ant mass is complementary to e.g. exclusive p- 
meson production and may offer some advantages 
because it is likely to be less affected by higher 
twist effects. Second, a quantitative understand- 
ing of dijet electroproduction is imperative in order 
to address more ambitious and theoretically more 
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FIG. 1: Kinematics of hard diffractive dijet production 



challenging cases of hard dijet pro duction by in- 
cidentpions 0, El El El El El and real pho- 
tons 1 13 . Ilfij . Another very important extension is 
exclusive dijet production in pp collisions, where 
one can study many aspects relevant for exclusive 
diffractive Higgs production ^ n the future, 

there might also be an interesting opportunity to 
study spin asymmetries in dijet production at eR- 
HIC accelerator, see e.g. E3- 

In theory, exclusive dijet production is defined 
as the process in which the invariant energy M 2 
deposited in the two narrow angular regions (jet 
cones), M 2 = (gi +Q2) 2 , almost coincides with the 

(p + q — p') 2 of the 



total invariant mass M 2 is 



diffractively produced system, M 2 > (1 — e)Al 2 iff , 
with e <C 1 serving as an infrared cutoff. In ex- 
isting experiments, diffractive events are usually 
defined by the presence of large rapidity gaps in 
the hadronic final state, and main observation that 
trigged the interest to hard diffraction has been 
that the probability of large rapidity gaps is not ex- 
ponentially suppressed. Somewhat imprecisely, we 
will refer to dijet production under the condition 
of a large rapidity gap in the hadronic final state 
as inclusive production. There are arguments that 
inclusive diffraction production of dijets with large 



2 



invariant mass is dominated by processes like the 
one shown in Fig. \T\) with a gluon (gluons) emitted 
from the gluon ladder (pomeron) in the i-channel. 
The gluon in Fig. ^> is emitted preferably in the 
central rapidity region, which corresponds to the 
case M| iff > M 2 . 

The experimental distinction between exclusive 
and inclusive dijets proves to be an intricate prob- 
lem. One possibility to make such a separation 
would be to study the corresponding event shapes, 
for example, by imposing a suitable cutoff in the 
heavy jet mass. Another proposal 01 i s to identify 
the exclusive dijet final state by requiring that the 
jet transverse momenta are compensated to a high 
accuracy within the diffractive cone and making 
some additional cuts. This approach seems to work 
for the case of coherent dijet production from nu- 
clei by incident pions, and for photoproduction the 
corresponding experimental program is proposed 
for HERA [llj. 

The dijet electroproduction at high Q 2 offers it- 
self as the simplest process of this kind, in which 
QCD collinear factorization can be established and 
allows one to make well defined and stable predic- 
tions for the exclusive dijet longitudinal momen- 
tum fraction distributions as well as the depen- 
dence on the azimuthal angle between the jet and 
the lepton planes; these distributions can be used 
to test the separation of exclusive dijets from the 
inclusive sample. The purpose of this paper is to 
work out the necessary cross sections and present 
detailed calculations for HERA kinematics. A simi- 
lar program was suggested in Ref . Q where a differ- 
ent theoretical approach based on fc t -factorization 
was used. Our calculation is also similar to []| 
where hard exclusive meson pair production has 
been considered. A comparison of this earlier works 
with our results is done below in the text. 

Throughout this paper we will work in a refer- 
ence frame where the virtual photon and the pro- 
ton collide back-to-back. We will neglect proton 
mass whenever possible, p 2 = p' 2 — ► 0, and choose 
the (almost) light-like incident proton momentum 
p^ to be in "plus" direction. It is convenient to 
introduce another light-like vector in "minus" di- 
rection: 

?m = , ?' 2 = 0, (q'-p) = (q-p) , (3) 

so that e.g. the lepton momenta can be decom- 
posed in the two light-cone components ~ p, q' and 
the orthogonal plane: 

l » = -% + -Pn + 1 ±» > 

V V 



1-2/ 

y 



(4) 



Note that l\ is defined as the square of the trans- 
verse plane vector, i.e. with opposite sign com- 
pared to the square of the corresponding four- 
vector. 

Further, let W be the invariant cm. energy of 
the virtual photon-proton scattering subprocess 

i*(q)p(p) -> q{qi)q{q2)p{p') , (5) 

i.e. s 7 . p = (q + p) 2 = W 2 . We define 

& = P'-P, P = V A^~, t = A 2 , 



(q-A) 2 = M 2 , x Bj = 



Q 2 



W 2 + Q 2 



(6) 



To our approximation M 2 = (qi +q2) 2 is an invari- 
ant mass of the diffractively produced system. 



II. KINEMATICS OF EXCLUSIVE DIJET 
PRODUCTION 

We introduce two light-like vectors 



n, = n 2 _ = , n + rt_ = 1 



(7) 



in such a way that 

p=(l + 0Wn+ =p+, 

* Q 2 + w2 , to, 

q= 2W(lW) n - = q -' (8) 

where £ is the usual asymmetry parameter [(| 
which defines, in the scaling limit, the plus com- 
ponent of the momentum transfer 



p+ -p+ 



(9) 



P++P+ 

Any four- vector a M is decomposed as 

a' 1 = a + n1+a~n>^+a^ , a 2 = 2 a+ a~ -a\ . (10) 
Then, in particular 
Q 2 + W 2 



p'={l-i)Wn + + 



2(1 ~0 W 



n_ + A ± . (11) 



We consider the case when jet transverse momenta 
are compensated, (qi + q 2 )±_ = -A ± = so that 
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the last two terms in the second equation in Ijllfl 
can be dropped. The momenta of the dijets can be 
written as 



Q 2 + W 2 
2W(l + 



zQ 2 

Q 2 + W 2 
2W(l + 

2 



zQ 2 



-W(l+Ox Bj n+-q ± , (12) 



where m is the quark mass and z is the relative 
longitudinal "minus" momentum fraction carried 
by the quark jet. Here and below we use a short- 
hand notation 



z . 



The dijet invariant mass equals 



( qi+q2 f = M 2 = q - +m 



zz 



(13) 



and the asymmetry parameter 

2£ _ Q 2 + M 2 
TTe ~ Q 2 + W 2 ' 

For future convenience wc introduce the notation 



(14) 



and the parameter 



9± 



m 2 + zzQ 2 . 



Q 2 



M 2 + Q 2 



(15) 



(16) 



which for light quark jets m — > coincides with the 
conventional /3-parameter used in the description of 
diffractive deep inelastic scattering, see e.g. 

We will calculate the distributions in the angle 
between the electron scattering and the jets planes, 
cf. 0, ^3 . In order to define this angle we intro- 
duce two transverse unit vectors 



e x =(0, 1,0,0) 
e y =(0, 0,1,0) 



(17) 



in such a way that the incident lepton transverse 
momentum Q is in x-direction, lj_ = |Zj_|e x , and 
the quark jet transverse momentum equals to q± = 
I <?_L | (e x cos <p + e y sin </>) , see Fig. |21 




FIG. 2: The azimuthal angle <f> is defined as the an- 
gle between the two planes: one formed by the lepton 
momenta I, I' and the other one by the jet momenta 

?1,92. 



The deep inelastic differential cross section is 
written as 

da=(2ir) 4 6 4 (l + p - I' - q x - q 2 - p') 

\M\ 2 d 3 l'd 3 q 1 d 3 q 2 d 3 p' 



4(1 -p) (2Tr) 12 2l' 2q lfi 2 q2 ,o2 P ' 



(18) 



The amplitude M can be expressed in terms of 
the amplitude of the photon-nucleon scattering. In 
the Feynman gauge we have 



M 



(19) 



and using gauge invariance can replace the ten- 
sor by the sum of projections on the different pos- 
sible polarization vectors of the virtual photon: 
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[At/ 



L 2^1 e T,\ e T 



T.X 



(20) 



In what follows we use the longitudinal polariza- 
tion vector 



o 2 + w 2 w 

n- + ^(l + 0x Bj n+,(2l) 



2WQ(l + 

and define the transverse polarization vectors as 
ct=Q ■ n + + • n_ + e± , 



e T,: 



(22) 



Accordingly, we define the photon subprocess am- 
plitudes for different polarizations as 



A L = M£ e L , M , At = e T ,» (23) 



and 



A) 



A\ = M^ey tll (24) 
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It is not difficult to show that 



and the cross section of interest takes the form 



d 3 l' 1 

w = - 4 dydQc 



(25) 



(27T) 



(1(7 



dydQ 2 d(j)i> nQ 2 y 



A-Ay + y 2 



,2 , V 



\A T \ T ^ |«T. T 



|^| 2 + (1 - y)|^| 2 + (2 - y)yi^i?e (A* T A* L ) 



dzd(f>dqj 



dA±d(f>p> 



2 9 n 4 zz(W 2 + Q 2 )(W 2 - M 2 ) 2tt 

I 



(26) 
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FIG. 3: Gluon contribution to the hard exclusive dijet FIG. 4: Quark contribution to the hard exclusive dijet 



production 



production 



Here </y is the azimuthal angle for the outgoing nu- 
cleoli. The corresponding integration is trivial so 
that the factor d(j> p i /2tt can be replaced by unity 
for all practical purposes. In the numerical cal- 
culation described below we assume the behavior 
dafdA\ ~ exp[— bA 2 ^} with the universal slope 
b = b GeV 2 . 



III. CALCULATION OF THE 
AMPLITUDES 

In this work we calculate the necessary ampli- 
tudes to leading order in perturbation theory us- 
ing collinear factorization. The corresponding di- 
agrams are shown in Fig. [3] and Fig. 0] for the 
gluon and quark contributions, respectively. In 
both cases the addition of symmetric diagrams with 
the permutation of the quark and the antiquark is 
understood. 

In the framework of collinear factorization, the 
necessary hadronic input is parametrized in terms 
of generalized parton distributions (GPDs). GPDs 
are defined as matrix elements of light-ray quark 
and gluon operators sandwiched between nucleon 
(proton) states with different momenta @, : 

F«(x,Z,A 2 ) = 



H"(x^,A 2 )u(p') i-u(p) 



+ £i(x,tA 2 )u(p') 
F9(x,£,A 2 ) = 



ix(Pz) 



2tt 



2 Tti 7v 



u(p) 



(Pn_) 



\z=Xt%- 



2(Pn_) 



io aS3 n_ a A f} 
2 tojv 



u(p) 



(27) 



Here u{p) and u(p') are the nucleon spinors. In 
both cases insertion of the path-ordered gauge fac- 
tor between the field operators is implied. The mo- 
mentum fraction x, —1 < x < \, parametrizes par- 
ton momenta with respect to the symmetric mo- 
mentum P = (p + p')/2. In the forward limit, 
p' = p, the contributions proportional to the func- 
tions £ 9 (x,£, A 2 ) and £ 9 (x,£,,A 2 ) vanish, whereas 
the distributions H q {x,£,A 2 ) and H 9 (x,£, A 2 ) re- 
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duce to the ordinary quark and gluon densities: 
H q (x,0,0)=q(x) for x > 0, 
H q {x,0,0)=-q(-x) for a;<0; 
H 9 (x,0,0)=xg{x) for x > . (28) 



Note that the gluon GPD is an even function of x, 
W{x,^A 2 ) = W{-x,^A 2 ). 

The calculation is relatively straightforward so 
that we omit the details. The 7*iV scattering am- 
plitude for the longitudinal photon polarization can 
be written as 



^7l= IgiT^p "wO^MUl + 2C F r L ) , (29) 

where 

lS -f dxF9(x e a 2 ) f 2 ^ + 2 ^ 2 ^ - ^ 

-1 
1 

11= / dxF q (x,^,A 2 ) ( 2 f + 2 f ) . (30) 
7 \{x + (i-ie) {x-£ + te)J 

-l 

We denote [3 — 1 — /3, u{q\) 1 v{q2) are the quark spinors and stands for the colors of the outgoing 
quarks, C F = (N 2 — 1)/2N C where N c = 3 for QCD. For the transverse photon polarization we obtain 



27ra s V47ra em e g % VF r _ „ 
+u{q x ) itix. ^+v{q 2 ){2C F I q T 1 + + fi(<n) fi + v(q 2 ){2C F I q2 - zI a T )} , (31) 



with 



£(1-2/3) £(1-2/3) 4£/3 



4= J dxF 9 {x^,A 2 ) 



-l 
l 

-l 

l 



(a; + £ — ie) 2 (x — £ + ie) 2 (x + £ — ie)(a; — £ + ie) 



If= I dxF q (x, £, A 2 ) ( ? - 5 , 2 ^" 2 , - + ^ 



(x - £ + ie) /3(x + £- ie) /3(x - £(1 - 2/3) - ie) 



I»- fdxFHx £ A 2 ) f ^ - 2 ^ 2 ^ (32) 

T "i l j l/3(x-£ + J£ ) (x + £-*e) ^ + £(l-2/3) + ie )J • (6Z) 

-l 

In both cases the expressions are written for one quark flavor. The quark mass corrections have been 
omitted for quark contributions, cf. Fig. 0] and only have to be taken into account in the charm- or 
eventually beauty-quark production in the gluon contributions shown in Fig. [3J Note that for heavy 
quark jets our definition of the /3-parameter (|16|) differs somewhat from the conventional /3-parameter 
used in the description of diffraction scattering. 

Using the virtual photon amplitudes given in (|29|l and 1311) we calculate the dijet cross section, summed 
over helicities and color of the produced (qq) pair: 



da a em 
dydQ 2 nQ 2 y 

where 



l + (l-j,) ; 



duT — 2(1 — y) cos2(f> ddTT + (1 — v)d^L ~ (2 — y)y/l — y cos ifidaLT 



(33) 



G 




FIG. 5: Longitudinal momentum fraction distribution of the dijets, summed over all quark flavors q = u,d,s,c. 
The three curves correspond to different minimum dijet transverse momenta; from above to below: qo = 
1 .25, 1.5, 1.75 GeV. The solid and the dashed curves correspond to a model for the generalized parton distri- 
bution based on the leading-order CTEQ6L and MRST2001LO leading-order parton distributions, respectively. 
The results shown on panel a) correspond to the calculation with a cutoff Q 2 /(Q 2 + M 2 ) > 0.5, whereas the ones 
shown on panel b) are calculated without such a cutoff. In the latter case we take Q 2 > 10 GeV 2 . 




b) 




FIG. 6: Azimuthal angle distribution of the dijets, q± > 1.25 GeV, summed over all quark flavors q — u,d,s,c. 
Identification of the panels and curves is the same as in Fig. EI see a l so text. 



d(JT 

d(JTT=Re 



^\I 9 L \ 2 + \zI 9 T + 2C F I qi | 2 + \zl* - 2C F I% 
1± 



92 |2 



da . 



(zI 9 T + 2C F If){zI 9 T - 2C F tf?y 



da , 



O 2 

d<r L =4(zz) 2 ^ T \I 9 L + 2C F I q L \ 2 da, 



da LT =2{zz)-^- Re (If + 2C F I q L ){{\ - 2z)l 9 T + 2C F {If + I 9 ?))* 



da 



(34) 



and 

a em a 2 s e 2 q q\ dq^d/^dzdcj) 
d(T -lQ*N c {l~e) (<zi+M 2 ) 4 • (5J 

The same cross section in the k j_ -factorization 
approach of Ref. 0] can be recovered by neglect- 
ing the double-pole terms in /£, Ij,, in (|3U|) . 



taking into account the imaginary part only 

/f^7r-2(l-2/3)^(e,e,A 2 ), 

I*->-t7r- W(£,£,A 2 ), (36) 

and also neglecting the quark mass and quark con- 
tributions altogether. The result coincides with the 
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expression given in |4|, except for the sign in the 
last term in ~ cos <j>. Our sign agrees with an 
independent calculation in 18] in a kind of a two- 
gluon exchange model. (To be precise, our defini- 
tion of <j) differs in sign from the definition used in 

IV. NUMERICAL RESULTS 

The numerical results presented below are cal- 
culated using a model for the generalized parton 
distributions that is based on CTEQ6L 20] (solid 
curves) and MRST2001LO [HI (dashed curves) 
leading-order parton distributions at the scale 
4 GeV 2 . For our analysis we neglect the contribu- 
tions of £ q and £ 9 , see l|27jl. since their contribu- 
tions is expected to be minor in the HERA energy 
range. The H q and H 9 distributions are modeled 
using the parametrization in terms of the so-called 
double distributions using an ansatz proposed in 
[22| with profile functions chosen to be tt(x, y) — 
6y(l — x — y)(l — x)~ 3 for quarks and ir(x,y) — 
30y 2 (l — x — y) 2 (l — x)~ 5 for gluons. We use the 
two-loop running QCD coupling corresponding to 
the value a s (Mz) = 0.1165. For the scale of the 
running coupling we take m 2 + q\ + zzQ 2 for the 
c-quark, and q\ + zzQ 2 for the light quark pro- 
duction. The charm quark mass is taken to be 
m c = 1.25 GeV. We plot one-dimensional differen- 
tial cross sections obtained by integrating over 
the remaining variables. The distributions over the 
longitudinal momentum fraction of the dijets ("z- 
distributions" ) are shown for the integral over all 
azimuthal angles, < <f> < 27r, whereas for the az- 
imuthal angle distributions ( "^-distributions" ) we 
integrate in the range 0.1 < z < 0.9. We take inte- 
gration limits in the deep- inelastic y- variable 0.1 < 
y < 0.4 which roughly corresponds to the energy 
interval W = 100-^200 GeV. We also integrate over 
the dijet transverse momenta q\ > q$ with three 
choices for the cutoff: qo = 1.25, 1.5, 1.75 GeV, and 
over Q 2 in the range Q 2 = 10 500 GeV 2 . In 
addition, we introduce a cutoff on the invariant 
dijet mass M 2 < Q 2 alias for the diffractive DIS 
/3-parameter 

/3 DDIS = Q 2 /(Q 2 + M 2 ) >0.5 

(cf. Qltifl h which is supposed to facilitate the ex- 
traction of the exclusive diffractive dijet cross sec- 
tion experimentally. Indeed, in the region of large 
/3 DDIS radiation of an additional gluon (gluons) in 
the final state, as shown in Fig. lb, represents a ra- 
diative correction 0(a s ) to the main process, the 
qq pair production. At the same time, for small 



/3 s radiation of gluons is enhanced by large log- 
arithms of the energy and is numerically very im- 
portant. It leads to events which have a topology of 
inclusive diffractive dijet production. In this case 
a special experimental procedure is needed to iso- 
late the exclusive contribution. Last but not least, 
it is worthwhile to mention that all calculations in 
this work refer to the parton-level cross sections, 
hadronization effects are not taken into account. 
We also average over the quark and the antiquark 
jets, since their distinction is very difficult experi- 
mentally. With this averaging, the last term ~ glt 
in drops out. 

The longitudinal momentum fraction and the az- 
imuthal angle distributions of the dijets, summed 
over all quark flavors q = u, d, s, c, are shown in 
Fig. [5] and Fig. respectively. Note that the 
z-distributions are affected strongly by the cutoff 
/jDDis > 0.5 [23, while the ^-distributions remain 
qualitatively the same. This effect is easy to under- 
stand and is due to a strong kinematic suppression 
of small z — > and large z — ► 1 longitudinal mo- 
mentum regions that correspond to large masses of 
the diffractively produced system. The difference 
between the solid and the dashed curves is mainly 
in the absolute normalization and it arises because 
of the different small- a; behavior of the CTEQ6L 
and MRST2001LO gluon distributions. 

We expect that the accuracy of our calculation 
is mainly limited by the size of the next-to-leading 
order (NLO) corrections. E.g. for vector meson 
electroproduction the NLO corrections were found 
to be large [24|. The uncertainties involved in the 
modeling of generalized (off-forward) parton distri- 
butions are probably less important in the HERA 
energy range that we consider in this work. 

The same distributions are plotted in Fig. [7| and 
Fig. [S] for the contributions of c-quark jets sepa- 
rately. Note that the cutoff (3 DDIS > 0.5 leads to 
a dramatic reduction of the cross section in this 
case, so that the c-quark contribution to the distri- 
butions in Fig. and Fig. is rather small. 

As it can be expected, the dijet production is 
dominated by the contribution of transverse pho- 
ton polarization, dar in Eq. The contribu- 
tion of the longitudinal polarization, daL, is shown 
separately in Fig. for the case q± > 1.5 GeV. 
It is seen that the relative weight of the longitu- 
dinal contribution is effectively enhanced by the 
cutoff /3 DDIS > 0.5 (since dax is peaked at small 
and large z and is strongly suppressed in the end 
point regions by the cutoff). 

In addition, in Fig. llOl we show the y-distribution 
of the cross section, integrated over the longitudi- 
nal momentum fraction of the jets, and over the 
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FIG. 7: Same as in Fig.|5] but for charmed quark dijets only. 
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FIG. 8: Same as in Fig. [HI but for charmed quark dijets only. 



azimuthal angle. This plot represent, essentially, 
the energy dependence of the cross section of vir- 
tual photon proton scattering (JSJ. In the consid- 
ered kinematic range we find a steep rise <7^*p ^ 
^2-j0.24h-o.26 wn j c h j g typical for hard diffractive 
processes. 

Finally, let us compare our results with the ones 
obtained earlier in 3, 4] within the fci-factorization 
framework. Qualitatively, both collinear factoriza- 
tion and fet-factorization approaches lead to simi- 
lar predictions: in particular the same overall en- 
sealing, strong energy dependence of the virtual 
photon proton scattering and the prediction that 
jets prefer a direction perpendicular to the elec- 
tron plane. This is not surprising since both tech- 
niques coincide in the double logarithmic approxi- 
mation, see l|36l) . The direct comparison of numeri- 
cal predictions is difficult due to different cuts used 
and also different parameterizations for the input 
gluon parton distributions. It seems, however, that 
our cross sections tend to be roughly factor of two 
smaller than the ones reported in [j| • 

Beyond the double logarithmic approximation 
there are differences. In collinear factorization 
quark GPD contributes together with gluon GPD 



at the leading order. It turns out that the quark 
GPD contribution to the amplitude is quite signif- 
icant. However, we find (as a representative exam- 
ple, see Fig. lll|l that it interferes destructively with 
the gluon GPD contribution and, as the result, the 
contribution of quark GPD to the cross section can- 
cels to a large extent with the gluon-quark GPDs 
interference term. Remarkably, this cancellation 
does not affect much the shape of different distri- 
butions, and only results in a moderate increase of 
the cross section, as compared to the calculation 
with only the gluon GPD taken into account. One 
noticeable exception is the longitudinal contribu- 
tion to the cross section calculated with a /3 DDIS 
cutoff. As seen in Fig. El neglecting quark GPD 
contribution underestimates the longitudinal cross 
section by more than a factor two. 

Another difference and a distinct feature of the 
collinear factorization approach is the appearance 
of double poles in the coefficient functions for the 
gluon contributions to both the transverse lp?|) and 
the longitudinal (|30[1 amplitudes. Such double-pole 
terms have no counterparts in the k t factorization 
approach. In Fig. ^] we show absolute values of 
I j, and I 9 L as defined in Eqs. (|32|l and (j20J), respec- 
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FIG. 9: Longitudinal momentum fraction distribution of the dijets with transverse momentum q± > 1.5 GeV, 
summed over all quark flavors q = u, d, s, c (the two upper curves). The contribution of the longitudinal photon 
polarization is shown separately (the two lower curves). The calculation is made with (panel a)) and without 



(panel b)) a cutoff /3 



Q /(Q + M ) > 0.5. The solid and dashed curves correspond to the CTEQ6L and 



MRST2001LO parametrizations, respectively. 
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FIG. 10: The y-distribution of the dijets, yda/dy, 
summed over all quark flavors q = u,d, s, c. The cal- 
culation is made with a cutoff Q 2 /(Q 2 + M 2 ) > 0.5. 
Identification of the curves is the same as in Fig. Qj] see 
also text. 



lively, calculated for a typical value of the asymme- 
try parameter £ = 0.001 © and for different values 
of the /3-parameter <|16|) . The double pole contri- 
butions (shown by dashed curves) are important in 
the regions where the main single-pole terms van- 
ish: at small f3 for the transverse amplitude and at 
f3 r*> 0.5 for the longitudinal amplitude. 

It is known that at large energies (small x) the 
amplitudes are predominantly imaginary. Also in 
our case the imaginary part of the amplitude dom- 
inates the dijet cross section. We find however, 
see Fig. that the contribution of the real part 
is quite sizable close to the end points, z — > or 
z — ► 1, and especially in the case of the calculation 
without a /3 DDIS -cutoff. 



FIG. 11: The longitudinal momentum fraction distri- 
bution of the dijets with transverse momentum q± > 
1.25 GeV summed over light quark flavors q — u,d,s. 
The contributions of the gluon and the quark GPDs 
are shown by the dashed and the dash-dotted curves, 
respectively. The sum of all contributions including 
the quark-gluon interference terms is shown by the 
solid curve. The calculation is made with a cutoff 
(3 Dms > 0.5 and CTEQ6L parton distributions. 

V. SUMMARY AND CONCLUSIONS 

We have presented a detailed analysis of exclu- 
sive diffractive dijet production with large trans- 
verse momenta in the framework of QCD collinear 
factorization. The calculation is done in leading 
order in the strong coupling. We derived the ex- 
pressions for the amplitudes for the qq pair produc- 
tion by the virtual photon both for the transverse 
and the longitudinal polarizations and used these 
results for an extensive numerical study of the dif- 
ferential cross section for HERA kinematics. We 
have checked that our results for the amplitude 
are equivalent in the massless quark limit to the 
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FIG. 12: Longitudinal momentum fraction distribution 
of the dijets with transverse momentum q± > 1.5 GeV, 
summed over all quark flavors q = u,d,s,c. The full 
calculation (solid curves) is compared with the result 
that includes the contribution of the gluon GPD only 
(dashed curves). The two lower curves show the lon- 
gitudinal contribution separately. The calculation is 
made with a cutoff f3 Dms > 0.5 and uses CTEQ6L 
parton distributions. 

ones obtained in [j| (where no separation into the 
transverse and longitudinal contributions is made 
and different notation used for the GPDs). In the 
double logarithmic limit our result agrees (except 
for the sign of ~ cos0 term in (|33[l ) with the one 
obtained in the kt factorization approach 4j. 

Experimentally, main challenge in the study of 
hard dijet production is the necessity to have a 
clean separation between the exclusive and inclu- 
sive channels. Since topology of the event is dif- 
ferent in these two cases, such a separation should 
be possible to achieve using appropriate cuts. The 
most practical possibility at present is probably to 
limit the study of diffractive dijet electroproduc- 
tion to the kinematic region of large /3 D , say 
/3 DDIS > 0.5, where the exclusive qq production 
represents the main contribution and radiation of 
an additional gluon (gluons) in the final state is 



suppressed. Our estimates indicate that in the re- 
gion /3 DDIS > 0.5 the cross section remains suf- 
ficiently large and we hope that such an analysis 
can be done at HERA. As first noted in Q , the az- 
imuthal angle distribution of the dijets can be used 
to check the separation of the exclusive sample: ex- 
clusive jets prefer a direction perpendicular to the 
electron scattering plane whereas in the inclusive 
case the distribution is peaked in this plane. We 
find that the azimuthal angle distribution is sta- 
ble to various cuts and is not very sensitive to the 
input GPDs, so it can indeed be used as a useful 
trigger. 

Though quark GPDs contribute significantly to 
the amplitude, we observe a large cancellation be- 
tween the square of the quark contribution and the 
gluon-quark GPDs interference term in the cross 
section. As the result, both the magnitude of the 
cross section and the shape of different distribu- 
tions appear to be not very sensitive to the pres- 
ence of quark GPD contributions (in the studied 
energy range). This finding further strengthens the 
existing expectations (see e.g. @) that exclusive 
diffractive dijet production may offer an interest- 
ing possibility to constrain the generalized gluon 
parton distribution at small Xb- 
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FIG. 14: Longitudinal momentum fraction distribution of the dijets, summed over the light quark flavors q — 
it, d, s. The three pair of curves correspond to different minimum dijet transverse momenta, from above to below: 
go = 1.25, 1.5, 1.75 GeV. The solid and the dashed curves correspond to calculations when both the imaginary 
and the real part of the amplitude and when only the imaginary part of the amplitude is taken into account, 
respectively. The results shown on panel a) correspond to the calculation with a cutoff /3 DDIS > 0.5, whereas the 
ones shown on panel b) are calculated without such a cutoff. In the latter case we take Q 2 > 10 GeV 2 . 
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